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L.G. Houpert and B.J. Hamrock 

E lastohydrod.vnamic Lubrication Calculations Used as a Tool to Study Scuffing 


A new Reynolds equation is developed that takes into account the nonlinear viscous behavior of the 
fluid. The new Reynolds equation considers the nonlinear viscous fluid model of Eyring, the equi- 
librium equation, the constant mass flow, and the kinematic boundary condition. The new Reynolds 
equation and the °lasticity equation are solved simultaneously by using a system approach and a 
Newton-Raphson technique. Comparisons are made with results obtained from the classical Reynolds 
equation. The effects of sliding speed and introducing a bump or a groove within the conjunction 
are studied. Results are shown for both moderate and heavy loads. 


1 INTRODUCTION 

Scuffing is a problem encountered mainly in gears 
and rolling-element bearings. In rolling-element 
bearings, for example, scuffing is found in the 
acceleration zone of large bearings that are 
heavily loaded. For these circumstances, inertia 
effects probably lead to high sliding speeds. 

The authors feel that scuffing is related to 

(1) Film thickness failure of the macro- 
contacts, as studied by Dyson (1976), 
end the microcontacts 

(2) High pressure and shear stresses act- 
ing on the surface 

(3) High surface temperatures found when 
sliding speeds are high 

Film thickness failures lead to metal-to-metal 
contact; high surface stresses lead to rolling 
bearing fatigue, as studied by Ioan.’ides an 
Harris (1985); and high surface temperatures m 0 y 
cause local melting of the surfaces. 

Since the surface temperature can be calcu- 
lated if the kinematic and surface stresses are 
known, important parameters to control when 
studying scuffing are film thickness and surface 
stresses. These parameters can be calculated by 
using elastohydrodynamic lubrication (EHL) ana- 
lysis. In the present paper, a new method is 
describee to study macro- and micro-EHL without 
any restriction on the load applied to the con- 
tact. By macro-EHL, we mean the lubricant film 
thickness developed in the inlet zone of the EHL 
conjunction. Micro-EHL may occur below the 
asperities and is due to the squeeze or sliding 
speed effects of the bearing as described by 
Cheng (1983). Macro-EHL is first studied by 
using the classical Reynolds equation (linear 
viscous fluid model), and the results are shown 
to be accurate for any applied loads. But under 
the severe conditions mentioned above (large 
pressure and high sliding speeds), the lubricant 
behavior can no longer be considered to be linear 
viscous, and a nonlinear viscous model is devel- 
oped. With this fluid model, a modified or new 
Reynolds equation is derived. The new Reynolds 
equation and the elasticity equations are solved 
simultaneously by using a system approach and a 
Newton-Raphson technique. Film thickness, pres- 
sure, and shear stress can be obtained without 
any load restriction. The method described 


herein is therefore a very powerful tool that 
can be used in studying scuffing. To illustrate 
the power of the new approach, stress concentra- 
tions near a bump and a groove are studied. 

1.1 Symbols 

A dimensionless sliding speed, 

(ub - u a )/u a 

b half Hertzian w 4 dth, R^stiih, m 

b' half-width of singularity, m 

Cj constant, 2E'W/-n 0>re * 

C? constant, iEU/8Wr 0>re f 

Cj weight factors used to integrate P 

E modulus of elasticity, N/m ? 

E‘ effective elastic modulus, 

2(11 . 

♦ (1 - ''b ) /E b )] _1 . N / m 2 

G dimensionless materials parameter, aE’ 

H dimensionless film thickness, 

hR/b 2 . ih/8RW 

H e dimensionless filir. thickness when dP/dX • 0 

Hg dimensionless constant used in calculation 
of H 

h film thickness, m 

1,J nodes 

K constant, 3» Z U/4W 2 

N number of nodes use'* in linear system 
P dimensionless pressure, p/p H 

p pressure, N/m 2 
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p maximum Hertzian pressure, E'b/4R, N/m 2 
H 

R equivalent radius, (l/r a ♦ l/rfo) - *, m 

r curvature radius, m 

S nu a /hr 0 ■ C 2 n/Hf 0 (1 + A/2) 

T (h/ To ) dp/dx . Ci(H/f„)(dP/dX) 

U dimensionless speed parameter, nou/E'R 

u velocity of fluid 

u mean velocity, (u a ♦ ) /2 

W dimensionless load, w/E'R 

w applied load per unit length, N/m 

X dimensionless abscissa, xR/b 2 

x abscissa along rolling direction, m 

x p position of singularity 

Z dimensionless coordinate, z/h 

z coordinate In direction of film thickness, m 

z m height of singularity, m 

a pressure-viscosity coefficient of fluid, 

nr/N 

t strain rate, s - * 

n viscosity of fluid, N s/m 2 

W dimensionless viscosity, n/no 

no viscosity at atmospheric pressure, N s/m 2 

w R traction coefficient due to rolling speed 

u S traction coefficient due to sliding speed 

v Poisson's ratio 

p density of lubricant, kg/m 3 

T dimensionless density, p/pq 

p e density when dp/dx > 0, kg/m 3 

P0 density at atmospheric pressure, kg/m 3 

t shear stress, N/m 2 

t o critical shear stress, N/m 2 

T 0 dimensionless critical shear stress, 
T o/ T o,ref 

2 CLASSICAL REYNOLDS SOLUTION 

From Houpert and Hamrock (1985), the Integrated, 
dimensionless form of the classical Reynolds 
equation may be written as 



where 



Conventional so', itlon of equation (1) required 
large computer run time, and the approach failed 
at moderate and high loads. Okamura (1982) used 
a system approach In solving equation (1). 

Houpert and Hamrock (1985) improved on Okamura' s 
approach by making the following changes: 

(1) Used a nonuniform mesh to Improve the 
accuracy of the calculations In the 
pressure spike region and In the Inlet 

(2) Used a more accurate way of calculat- 
ing the elasticity 

(3) Solved for pressure difference rather 
than actual pressures 

(4) Used Roelands' pressure-viscosity model 
By incorporating these changes, Houpert and 
Hamrock (1985) were able t w get accurate and 
fast solutions without any restriction on load. 
Figure 1 Illustrates how fast the new approach 
is. It gives the pressure and film profiles at 
Iterations 0, 1, and 14. Note that In the first 
Iteration a pressure spike was formed that Is 
close to the final converged pressure spike. In 
general. It took about 15 Iterations to obtain a 
converged solution, or about 2 min of CPU time 
on an IBM 370 computer with a mesh of 181 nodes. 
Compare this with the 100 min normally taken by 
the Hamrock and Jacobson (1984) approach. 

Figure 2 shows the pressure profile and film 
thickness ratio for five dimensionless loads 
varying over two orders of magnitude. This 
dimensionless load range corresponds to a maxi- 
mum Hertzian pressure of 0.4 to 4.8 GPa, which 
is well within the range that rolling-element 
bearings and gears experience. As the loau 
increases (Fig. 2(a)), the pressure spike becomes 
smaller and moves toward the outlet. Further- 
more, as the load increases, the inlet meniscus 
moves toward the abscissa X ■ -1. The nip film 
thlcxness width and length (Fig. 2(b)) both 
decrease as the load Increases. 


3 NEW REYNOLDS EQUATION 

The classical Reynolds equation (1) has been 
obtained by assuming the lubricant behavior to 
be linear viscous. This Is only true at low 
pressures, or high temperatures, and when the 
shear rates are small. A more appropriate lubri- 
cant rheological model Is the Eyrlng (1936) non- 
linear viscous model used by Hirst and Moore 
(1974), Johnson and Greenwood (1980), and Berthe 
et al. (1978 and 1979). Elastic viscous models 
have also been used by Johnson and Tevaarwerk 
(1977), Bair and Miner (1978), Berthe et al. 
(1979), Houpert et al. (1981), and Houpert 
(1985a), but elastic effects have been shown to 
be often negligible as proposed by Houpert 
(1985b). 

The Eyrlng (1936) model will therefore be 
used to develop a new Reynolds equation. The 
new Reynolds equation will assume *n Isothermal 
K ehav1o> and neglect transient and squeeze-film 
eti However, the new rheological model. 
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equilibrium equation, boundary kinematic condi- 
tions, and constant mass flom conditions have 
been taken Into account. 

From equilibrium, we can write that 


*■ 


31 


(3) 


Integrating gives 


Ir. T „ dx T„ 


(4) 


where A} Is a constant to be defined later. 
Erylng's (1936) nonlinear viscous fluid model 
can be expressed as 


du 

Y m Tz 


■ T >'" B fe) < 5 > 

Substituting equation (4) Into equation (5) gives 


( 6 ) 


Integrating gives 




♦ A, 


where A? Is an Integration constant. The 
boundary conditions are 


(7) 


(1) z - 0, u - u s 

(11) z - h, u - u h 


By making use of these boundary conditions, the 
constants A} and A 2 are defined and we 
find that 

H- . 1 . ^ [co,* (zr • - c.,» (gj 


where 


*1 -1 
— - slnh 

0 


u K - u. _ 

. b a , z 

A - — uT - ' Z ‘ T>’ 


AST 


2 slnh 




T 

7 


(9) 


nU a 

“ tTF “ 


FI 


( 10 ) 


T h dp r H dP 
T ” 7“ 37 " C 1 = 37 
0 T o 

Making use of equations (9) and (10), we can 
rewrite equations (4) and (6) as 


T - > slnh* 
T o 


l[_AST_^ 
I 2 slnh J 

■ °' 5) 

i 


♦ T(Z - 0.5) (11) 


♦ slnh 


AST 


2 slnh 

T* 

7 


(12) 


To calculate the mass flow. It Is convenient 
to express the velocity distribution as 

s-.i.n.f,. .!■*[} ;z-i)] 


x slnh j-^ * slnh** 


AST 


2 slnh ^ 

5 


(13) 


Note from equation (13) that when dp/dx * 0, 
then T* 0, slnh (T/2)/(T/2) * 1, and 
slnh [T/2(Z - l)]/(T/2) * (Z - 1). Therefore 


. 1 ♦ AZ when T * 0 
u a 


The mass flow when dp/dx - 0 can be 
written as 


j[ V. (S;) « ■ *A | I* * «> « 

• p e h e (* * f) 

The condition of constant mass flow at any 
location can therefore be written as 

•A (> * 7 ) - *" (ij) « 


(14) 


(15) 


(8) Making use of equation (13) gives 




1 ♦ 


AST 


2 slnh 


[f slnh (J) - cosh (J) 


1/2 


(16) 


Equation (16) Is the new Reynolds equation. But 
In order to make a comparison with the linear 
viscous equation, equation (16) must be rewritten 
as 


j3 dP 


37-7 





' AST ' 

2 

1 ♦ 

t 


slnh ^ 



117) 
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where 



Comparing equation (1) with equation (17) Indi- 
cates quite clearly the contribution due to the 
norllnear fluid model . 

Table I shows the difference In the equa- 
tions developed for the nonlinear and linear 
viscous fluid models. Note that If t«* «, 
then T * 0, e * 0, and S * 0. For this situa- 
tion the nonlinear viscous equation reduces to 
the linear viscous equation. Therefore, the 
new Reynolds equation to be solved at each mode 
In dimensionless form can be expressed as 


appendix. The unknowns of the system are aPj, 
A(5 e H e ), and aHq; they must be less than or 
equal to, respectively, 1/1000 of Pj, 5 t H», 
and Hq. More details on handling tne exit 
boundary and. In general, how the system approach 
Is used are covered In Houpert and Hamrock (1985) 
and will not be repeated here. 

In general, 10 to 15 iterations are required 
for convergence of the results and only 2 min of 
CPU time are necessary on the IBM 370. This 
approach Is therefore very fast, powerful, and 
relatively easy to program. As a result of the 
calculations, the film thickness Hi and pres- 
sure Pi are known at every node. Further- 
more, the shear stress t, the shear rate 
»u/»z, and the velocity u are known for values 
across the film. 



4 SYSTEM SOLUTION 

The basic equation to solve at each node I Is 
the new Reynolds equation given In equation (19) 
or fi ■ 0, where f Is a function of p e H e , 

Hq, «od the pressures f>2 t0 P N* The expres- 
sion for f also contains the dimensionless film 
thickness H, which is expressed as 

|| 

H i- H o + r- + £ d ij p j "f? 9(x) (20) 

where g(x) defines the microgeometry. The 
numerical approach used to solve the system of 
equations fi - 0 is Identical to the one used 
by Houpert and Hamrock (1985). A Newton-Raphson 
technique is used and leads to the solution of a 
linear system of equations where a Jacobian 
matrix has to be calculated as a function of the 
Jacobian elements »fi/»(5 e H e ), afj/aPj, and 
afi/»HQ. A linear system of N ♦ 1 equations 
Is solved 


” ,f l 

‘ f l 

^ • 

* f l 

• ^ 

•'7 

w o 


/*<W 


f-M 

,f 2 

•fW 


,f 2 

• ^ 

• f Z 


r | 


1 / 

,f N 

* f N 

Tir • 

• f N 

•• 

"5 


* “V ^ 

• 

1 

" f N 

0 

C Z • 

• C N 

0 


*<v 1 


»w / 


( 21 ) 

where C< are the Integration coefficients 
developed In Houpert and Hamrock (1985) and aW 
is a parameter representative of the constant 
load. The Jacobian elements are given In the 


5 DISCI SSI ON AND RESULTS 

Maintaining a fluid film of adequate magni- 
tude is an essential feature for the correct 
operation of lubricated machine elements. 

Results presented In this section make use of 
the nonlinear viscous rheological model of the 
lubricant and apply surface Irregularities such 
as Lumps and grooves to Illustrate the possibil- 
ities of the new approach. 

He first determined whether there was any 
difference between the film thickness obtained 
from the new Reynolds equation and that obtained 
from the classical Reynolds formulation. From 
disk machine experiments and traction force 
calculations, to has often been found to 
Increase linearly with pressure. But at high 
pressure, a limiting value t Is fourd 
as proposed by Houpert (1985b)’and Ten Napel 
et al. (1985). A possible relationship can be 
written as 


Jo_ . T o,ref T o . j . exp / _«E_\ ♦ , 

T 0,l T 0,l \ T o,i/ 

where a and i 0 t are defined for each 
lubricant as a function of the temperature. The 
term c Is very small ( 10 — 3 ) and Is Introduced 
to avoid numerical problems when calculating 
1 /t 0 as p * 0. 

Numerical runs have been performed with 
a w 2.6x10" and t q ref “ *o t * 1*^ JO Pa 
that show some small’numerical Instabilities at 
the outlet of the contact while the film thick- 
ness is essentially a function of the Inlet as 
clearly indicated by Houpert and Hamrock (1985) 
in their pressure spike analysis. The results 
obtained Indicate clearly that the film thick- 
nesses for the two forms of the Reynolds equation 
were Identical for.the conditions simulated here 
(l.e., U - lxlO" 11 , G . 5007, and M . 2xl0" 5 ). 
The reason for this Is that the shear stress was 
smaller than t 0 In the Inlet zone. Since 
the film thickness is established by the hydro- 
dynamic action In the inlet region and this was 
not altered by the new form of the Reynolds 
equation, the film shape did not change much. 
Because of this result and to avoid numerical 
Instabilities, t 0 will be assumed to be 
constant (t 0 - t 0 . ; t q - 1) In what 
follows But it should be mentioned that, In 
general, accurate film thickness calculations 
require the knowledge of t 0 at very small 
pressures although t 0 has always been 
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determined at large pressure for traction forces 
calculations. 

The effect of dimensionless sliding speed 
A on pressure and film shape through the lubrl- 
ca‘,ng conjunction Is shown In Fig. 3. The 
sliding speed had little effect on film thick- 
ness but a significant effect on the definition 
of the pressure spike. As the sliding speed 
increased, the pressure spike diminished. The 
main reason for the change In pressure spike 
with sliding speed was that near the pressure 
spike the shear stres* was larger than t 0 . 

The rolling traction coefficient u R remained 
almost unchanged (u R ■ 4.4x10“*! while the 
sliding traction coefficient n S increased 
from 0 to 8.4xl0“ z when A Increased from 0 
to 10. 

Ha/lng shown some results of using the new 
Reynolds equation, the nett task was to prestudy 
mlcro-EHL. As a first step and example of appli- 
cation of thw tool presented herein, we were 
interested In what effect a bump or groove has 
on the film shape and stress concentration. The 
change of shape due to a bump or groove Is taken 
care of by g(x) 1" the film shape equa- 
tion (20). For simplicity, g(x) is a polyno- 
mial of degree 2 described as a function of the 
singularity half-width b'. position x p , and 
maximum height z m , which is positive for a 
groove and negative for a bump. 



In what follows, b' Is equal to 0.25x10“* m 
and z m to -0.1x10“° m for the bump and 
2x10"“ m for the groove. These values are to 
be compared with the equivalent roller radius 
R - 0.111125xl0~ 1 m and the Hertzian half-width 
b - 0.8x10"* m at moderate load (W - 2xl0“ 5 ). 

The central film thickness value h in a smooth 
situation Is 0.25x10“” m. A large sliding 
speed, A = 10, has been imposed. 

Figure 4 shows the film shape and pressure 
profile when a bump occurs In the Inlet region 
of the conjunction (x B - -O.lxlO -3 m « -1.25 b). 
The undeformed and deformed shapes of the bump 
are presented. This figure shows that the bump 
occurring In the Inlet region had quite an effect 
on the level of the film thickness but did not 
alter the film shape appreciably. The film 
thickness was less with the bump than without. 

The bump also moved the pressure spike toward 
the outlet and lowered the spike height. The 
coefficients u R and were slightly 
affected by the bump ( W R . 6x10"* and 
w s . 8.56xl0 -2 ). 

Figure 5 shows the effect of a bump In the 
center of the conjunction on the pressure and 
film shape. The bump was the same as that used 
In Fig. 4, but Its location was changed (x p - 0). 
When the bump was In the center, there was no 
difference in minimum film thickness. Recall 
from Fig. 4 that the bump did significantly 
affect film thickness when It was located in the 
Inlet region. Furthermore, Fig. 5 shows that 
the presence of the bump in the center of the 
contact had a profound effect on the pressure 
profile. Stress concentrations p/pn of the 
order of 1.5 were calculated. Also the shape of 
the bump changed considerably In going from the 


undeformed state to the deformed state. This 
was not found when the bump was located In the 
Inlet region (Fig. 4). The coefficients u R 
and u 5 , remained almost unchanged (u R - 
4.9x10 , m S - 8.5xl0- z ) . 

The effect of a groove Is shown In Figs. 6 
and 7. The depth of the groove had to be made 
20 times that of the bump used In Figs. 4 and 5 
In order to get a comparable effect on pressure 
and film shape. Therefore, a groove would have 
considerably less effect on pressure and film 
shape than a bump If the bump and groove each 
had the same depth and width. Figure 6 shows 
the groove In the Inlet, and Fig. 7 shows the 
groove In the center of the contact. The obser- 
vations made for a bump In the Inlet and In the 
center of the contact can also be made about a 
groove, but substantial differences were noted 
concerning P R and u 5 . When the groove was 
In the Inlet, u R - 7.5x10"* and - 9*6x10-2 
as compared with u R - 3.98xl0“ 3 and u 5 - 7.89xl0“ z 
when the groove was In the center. Although 
temperature effects have been neglected In this 
study, the increase of v 5 when the groove 
was In the Inlet can be easily related to the 
experimental results of Cusano and Wedeveo 
(1983), who found similar Increases In u • 

The results presented thus far are for a 
moderate load. Figures 8 and 9 are for a heavy 
load (W - 3xl0“ 3 ). At this heavy load, 
b - 0.97xl0“ 3 m and h ■ 0.89x10““ m; 
b' - 1.94x10"! m (or 0.2 b) and 
z m = 1.48x10“!! m (or 16.6 h) for the groove 
and -1.43x10“' m (or -1.6 h) for the bump. 

The sliding speed A was fixed at 0.1. Figure 8 
shows the effect of a bump and Fig. 9 shows the 
effect of a groove, respectively, at positions 
Xp of 0.4 b and 0.2 b. The shapes of the bump 
and groove changed considerably once elastic 
deformation was considered. The pressure profile 
for the bump was very much Hertzian except at the 
corners of the bump. Stress concentrations were 
calculated for the groove only, but it should be 
r< -lembered that the groove was 10 times deeper 
than the bump. The coefficient U R was found 
to be almost nil (u R . 6.6x10“’) while 
was found in the two cases to equal 11.4; 10“ z . 

This large value was due to the isothermal 
assumption. 

6 CONCLUSIONS 

A new Reynolds equation was devf loped that takes 
Into account the nonlinear viscous behavior of 
the fluid. The new Reynolds equation and the 
elasticity equation were solved simultaneously 
by using a system approach and a Newton-Raphson 
technique. The film thickness was Identical to 
that calculated by using the classical Reynolds 
equation (assuming a linear viscous fluid model). 

The pressure spike decreased as the sliding speed 
increased. A bump or groove In ‘ e inlet had a 
more significant effect than If It was located 
at the center of the contact. The effect of a 
bump on pressure and film thickness was more 
significant than that of a groove. At heavy 
loads the impressive result was that, when elas- 
tic deformation was considered, the bump or 
groove essentially disappeared. These prelimi- 
nary results show that mlcro-EHL can be studied 
easily by using this new powerful tool. An 
exhaustive study can now be performed to control 
the effect of load, sliding speed, and roughness 
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wavelength and height on the stress field. Tem- 
perature Increases In the film and on the sur- 
faces should also be Included to better 
understand scuffing. 

APPENDIX 

Derivation of Jacobian Factors 

Calculation of df^/dPj. - Roelands' 
v 1 scoslty 1$ expressed as 

"n^ - exp {(In *0 ♦ 9.67) 

x [-1 ♦ (l ♦ 5.1x10”® 

where Z' Is the Roelands pressure parameter, 
equal to 0.69 In this study. By writing now 

dnj/dPj . A, we have 


*o « 


T 0,l \ T 04 / 


1 J 


Using the dimensionless variables we have 
dP 


T - C 


3x 


H 


1 - 


1 + 7 *o 


which leads to 


dT 1 ; c [( PH + "S) T o” T oS H ] 

■m- - T - C, * ■* 


n - 5.1x10"® p H [ln hq ♦ 9.67] 

*7 1 


-2 

T n 



r 

I 

* [l ♦ 5.1x10”® p^J n 1 k 1d 

dS i ; C 2 

nHT p 

- " « - 1 

where k^ Is the Kronecker symbol 

dP j l ♦ J 

CM 

£ 
i 


‘Ij " * " ' J 

k 1j - 0 If 1 a j 


The density o Is written as 


We can now define a series of new variables 
AST 


B 


2 sinh 


19' 


d(B.) . 

~3FT m B 1 


r j 


■ 1 ♦ 


0.6x10”® p H P 1 
1 ♦ 1.7x10'® p^ 


R A 

B 1 “ 7 


which leads to 
dp^ 

W~ ia 


0.6x10"® p u 


7 k 1j 


i (l ♦ 1.7x10'® Ph P^ 

In the appendix of Houpert and Hamrock (1985), 
d ( dP/dX ) /dPj has been expressed and can be 
written for simplicity as 


(ST + ST) sinh - T cosh Q- 

)" 

[ 1 

sinh (J)]* 


2 \ 1/2 

d(B ? ) . 

B 1 B 1 

+ B l) ; 

dPj ‘ 8 2 ; B 2 ‘ 

"V 

sinh 
B 3 — r 

^ - cosh (J) ♦ y? 



,/dP\ 

V). 


By definition of Dy we have 
dH. 

37~ ‘ H " D 1j 

For simplicity, the subscript 1 will be 
omitted when possible in what follows. 
From the definition of T 0 we have 


dB, 


i “ _ (II ? 7 — . ( i) . * ,» (}) . 


r - 12H „ . 
B 4 " ~3P~ B 3* 
HT 


«*(b 4 ) 
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TABLE I. - IMPORTANT EQUATIONS FOR NONLINEAR VISCOUS AND LINEAR VISCOUS MODELS 






















DIMENSIONLESS PRESSURE, P - plp H j AND FILM SHAPE, H/2H 




X-COORDINATE, x/b 

Fig. 1 Pressure and film profiles at iterations 0, 

1, and 14. Barus 1 pressure-viscosity formula. 
Dimensionless load, speed, apd materials param- 
eters held fixed at W • 2xl0" 5 , U * lxlO"U, and 
G ■ 3007. (From Uoupert and Hamrock (1985). ) 
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Fig. 2 Dimensionless pressure and film profiles 
for varying dimensionless loads. Roelands' pres- 
sure-viscosity formula. Dimensionless speed and 
materials parameters held fixed at U - lxlO' 11 
and G - 3007. {From Houpert and Hamrock (1983). ) 


DIMENSIONLESS 



X -COORDINATE, x/b 


Fig. 3 Dimensionless pressure and him profiles 
for varying dimensionless sliding speeds. 
Dimensionless load, speed, and materials 
parameters held fixed at W ■ 2X10*5, 

U - 1x10*11, and G - 5007. 



DIMENSIONLESS PRESSURE, P - p/p H ; 
AND FILM SHAPE, H/21HJ 



Fig. 4 Dimensionless pressure and film profiles 
without and with a bump in inlet region. Di- 
mensionless loaJ, speed, materials, and sliding 
speed parameters held fixed at W = 2xl(f 5 , 

U - lxlO" 11 , G - 5007, and A ■ 10. Bump 
depth, 0.1xl0“ 6 m; bump width, 0.5xl0"4m. 



DIMENSIONLESS PRESSURE, P - pfp H ; AND FILM SHAPE, H/2(H e ) 



X -COORDINATE, x/b 


Fig. 5 Dimensionless pressure and film profiles 
without and with a bump in center of contact 
Dimensionless load, speed, materials, and slid- 
ing speed parameters held fixed at W - 2x10^, 

U - lxlO-n, G - 5007, and A - 10. Bump 
depth, aixl0' 6 m ; bump width, 0.5xl0'4ni. 


DIMENSIONLESS FILM SHAPE, H/2(H e ) 



Fig. 6 Dimensionless pressure and film profiles 
without and with a groove in inlet region. Di- 
mensionless load, speed, materials, and sliding 
speed parameters held fixed at W ■ 2x10"^, 

U ■ lxKT 1 *, G - 5007, and A ■ 10. Groove depth, 
2xl0" 6 m; groove width, 0.5xl0‘ 4 m. 



DIMENSIONLESS FILM SHAPE, H/21HJ 


PROFILE WITH GROOVE 
PROFILE WITHOUT 


X-COORDINATE, x/b 

Fig. 7 Dimensionless pressure and film profiles 
without and with a qroove in center of contact 
Dimensionless load, speed, materials, and slid- 
ing speed parameters held fixed at W - 2xl0~ 5 , 

U ■ lxlfT^V G ■ 5007. and A - 10. Groove 


U ■ lxlO -11 , G - 5007, and A - 10. Groove 
depth, 2xl0‘ 6 m; groove width, 0.5xl0“ 4 m. 


DIMENSIONLESS PRESSURE, P 






DIMENSIONLESS PRESSURE: P * p/p H ; AND FILM SHAPE, H/2(H e > 


1.50 



X-COORDlNATE, x/b 

Fig. 8 Dimensionless pressure and film profiles with 
bump in heavily loaded contact Dimensionless load, 
speed, materials, and sliding speed parameters held 

fixed at W • 3xl0" 3 , U - lxlO' 11 , G - 5007, and 
A -0.1. Bump depth, 1.4xl0’ 7 mj bump width, 

3. 88xl0- 4 m. 


DIMENSIONLESS PRESSURE: P - p/Pns AND FILM SHAPE. H/2CH e ) 



Fig. 9 Dimensionless pressure and film profiles 
with groove in heavily loaded contact Dimension- 
less load, speed, materials, and sliding speed 
parameters held fixed at W - 3xl0' 3 , U - lxlO" 11 , 
G - 5007, and A - 0. 1. Groove depth, 1. 48xl0“ 6 m 
groove width, 3.88xl0" 4 m, 
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